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ABSTRACT 

Very high energy gamma-ray emission from the BL Lac object Markarian 421 has 
been detected over three observing seasons on 59 nights between April 1992 and June 
1994 with the Whipple 10-meter imaging Cherenkov telescope. During its initial detec- 
tion in 1992, its flux above 500 GeV was 1.6xl0~^^photons cm^^ s^^. Observations in 
1993 confirmed this level of emission. For observations made between December 1993 
and April 1994, its intensity was a factor of 2.2±0.5 lower. Observations on 14 and 15 
May, 1994 showed an increase over this quiescent level by a factor of ~10 (Kerrick et al. 
1995). This strong outburst suggests that 4 episodes of increased flux measurements on 
similar time scales in 1992 and 1994 may be attributed to somewhat weaker outbursts. 
The variability of the TeV gamma-ray emission from Markarian 421 stands in contrast 
to EGRET observations (Lin et al. 1994) which show no evidence for variability. 
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1 University of Michigan, Ann Arbor, MI 48109 
^University of Leeds, Leeds LS2 9JT, UK 

■^Whipple Observatory, Harvard-Smithsonian CfA, Amado, AZ 85645 

*Iowa State University, Ames, lA 50011 

^St. Patrick's College, Maynooth, Co. Kildare, Ireland 

^University College Dublin, Belfield, Dublin 4, Ireland 

^Purdue University, West Lafayette, IN 47907 

^Postal Address: Department of Physics and Geology, Northern Kentucky University, Highland Heights, KY 41099 
^Postal Address: Instituto de Astronomia y Fi'sica del Espacio, C.C. 67, Sue. 28, 1428 Buenos Aires, Argentina 



1 



1. Introduction 

The central engines of active galactic nuclei (AGN) 
are presumed to be massive black holes whose accre- 
tion powers radiation over a wide range of the elec- 
tromagnetic spectrum. More than a score of these 
have now been reported as sources of multi-hundred 
MeV/GeV photons (Fichtel et al. 1994), all of which 
are identified with the sub-class of AGN which are 
radio-loud and core dominated, and most of which 
are blazars. (See Dermer & Schlickeiser 1992, for a 
glossary of terms relevant to AGN). 

The AGN reported by EGRET vary in brightness 
by a decade and range in redshift from 0.031 to above 

2. The nearest such object is Markarian 421, a BL 
Lac object extensively observed at radio (Owen et 
al. 1978; Zhang & Baath 1990), UV/optical (Maza, 
Martin, & Angel 1978; Mufson et al. 1990), and X- 
ray frequencies (Mufson et al. 1990; Mushotzky et 
al. 1979; George, Warwick, & Bromage 1988). Its 
parent galaxy has been identified as a giant ellipti- 
cal (Ulrich et al. 1975; Mufson, Hutter, & Kondo 
1989) and it was the first BL Lac established as an 
X-ray source (Ricketts, Cooke, & Pounds 1979). Us- 
ing VLBI techniques, the radio fine structure of the 
source has been mapped with a resolution of 1 mas to 
show a core-jet-like structure which exhibits possible 
superluminal motion (Zhang & Bhaath 1990). Vari- 
ability with a time scale of hours has been observed 
in X-rays (Giommi et al. 1990) and at lower energies 
(Xie et al. 1988). Historically the source has exhib- 
ited low amplitude variability on the scale of months 
to years in the radio emission (Aller & AUer 1995) 
but so far there is no evidence for variability in the 
gamma-ray emission in the high energy regime (100 
MeV - 10 GeV, Lin et al. 1994). 

A comprehensive summary of recent theoretical 
work relevant to the blazars is given by von Mon- 
tigny et al. (1995). Although the ultimate origin of 
the blazar power is the central engine, high energy 
gamma-ray emission is likely to be beamed from a 
jet of highly relativistic particles. (If one assumes 
isotropy for the emission, enormous luminosities, in 
some cases more than 10^^ erg/s, would be required.) 
Furthermore any gamma rays at GeV-TeV energies 
that may be created in the inner regions of the source 
would be readily absorbed via the pair production 
process. Because of the likely jet geometry associated 
with the high energy emission process, time varia- 
tions observed at a single energy, if present, do not 



place strong constraints on the region of emission 
since there are many ways in which that variability 
may arise. However, correlated time scale variations 
at different energies (e.g. hard synchrotron and TeV) 
are generally useful and are the predicted consequence 
of models in which the gamma rays are produced by 
inverse Compton scattering. 

It has been pointed out by Gould and Schreder 
(1967) and recently reformulated by several authors 
(Stecker, De Jager, & Salamon 1992; Salamon, Steck- 
er & De Jager 1994; Dwek & Slavin 1994; MacMinn 
& Primack 1995) that the interaction of gamma rays 
with intergalactic infrared and optical photons pro- 
duced by galaxy formation may severely attenuate the 
TeV photon emission from cosmologically distant ob- 
jects and that evidence for such absorption (or lack 
thereof) may be used to constrain both the energy 
density of starlight and the distance of the sources. 
The dynamic energy range covered by the Whipple 
Observatory instrument is particularly sensitive to 
such effects for the assumed distance to Markarian 
421 (124 Mpc, Ho=75 km s~^ Mpc~^) and an upper 
limit on the IR density field from our observation is 
derived in a separate paper (Biller et al. 1995). 

2. Instrument and Observations 

The observations reported in this paper were car- 
ried out with the very high energy gamma-ray tele- 
scope (Cawley et al. 1990) at the Whipple Observa- 
tory on Mt. Hopkins in southern Arizona (31° 41' 
north, 110° 53' west, ah = 2300 m). This 10 meter 
diameter telescope images Cherenkov light from air 
showers on a hexagonal array of 109 fast photomulti- 
plier tubes (PMT's) in the focal plane. The inner 91 
2.86 cm diameter tubes of this high resolution camera 
have a 0.25° spacing. A surrounding outer ring of 18 
5.08 cm diameter PMT's was replaced in August 1993 
with a partial ring of 18 2.86 cm tubes. 

The effective energy threshold Eth, above which 
there is 100% detection efficiency over the effective 
collection area, was determined to be 0.4 TeV by sim- 
ulations for the 1988/89 observing season. For later 
periods, an adjustment due to changes in mirror re- 
flectivity and camera configuration was made by cal- 
ibrating with the measured cosmic ray rate R assum- 
ing a scaling behavior as i? oc E^^^ '^ for the energy 
threshold Eth- An energy threshold of 400 GeV cor- 
responds to an average count rate of 270 min"-'^. This 
scaling method can also be used to estimate gamma- 
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ray event rates at different energy thresholds for a 
source with a power-law spectrum and known 

spectral index a. In the context of this paper, we 
assume for the photon emission from the Crab Neb- 
ula an integral energy spectrum of the form E~^-^ as 
derived by Vacanti et al. (1992) 

A shower image is recorded when at least m out 
of the inner 91 tubes exceed the hardware threshold 
of n photoelectrons. With the standard trigger con- 
dition (m=2, n=50) a typical trigger rate of 4-5 Hz 
with i^t^ =400 GeV is obtained. All the observations 
were restricted to zenith angles less than 40° where, 
on the basis of our Crab Nebula observations, the 
rate variation with zenith angle is smaller than 10%. 
More detailed descriptions of the instrument and the 
Cherenkov light imaging technique can be found else- 
where (Cawley et al., 1990, Lewis 1990, Reynolds et 
al., 1993). The majority of the data were taken in 
a standard on/off tracking mode in which the source 
region and a background region are observed for equal 
times (usually 28 minutes). The location of the back- 
ground region differs from that of the source region 
only in right ascension; thus both regions have the 
same zenith angle coverage. This procedure elimi- 
nates systematic effects introduced by the small vari- 
ation in counting rate due to zenith angle changes. 

For a well established source such as Markarian 
421, the effective observing time can be doubled by 
observing in an alternative mode in which the tele- 
scope always follows the source ('tracking mode'). 
Under these conditions, shower selection procedures, 
described below, are used to establish an appropri- 
ate value for the rate of background showers. Table 1 
summarizes the varying detector parameters for the 
three observing periods between 1992 and 1994. 

Data were processed following the standard Whip- 
ple analysis procedures (Reynolds et al. 1993) in 
which the Cherenkov light images arc flatficldcd, 
cleaned and characterized by simple image moment 
parameters (Hillas 1985). The shape of the approx- 
imately elliptical shower image and its location rela- 
tive to the assumed source location in the image plane 
are specified by the parameters width, length and dis- 
tance. The fourth parameter, alpha, is defined to be 
the angle between the major axis of the shower image 
and a line from its centroid to the assumed source 
location in the image plane. For gamma-ray show- 
ers from a point-like source, alpha should be near 0° 
since the elliptical images point to the location of the 
source in the image plane. 



Information on the shape and the orientation of 
the images is used to remove more than 99% of back- 
ground cosmic-ray showers, while keeping a substan- 
tial fraction > 0.5) of possible gamma-ray showers. 
For this analysis we applied the same discrimina- 
tion values as we have for the selection of gamma- 
ray events from the Crab Nebula {0.073° <width< 
0.150°, 0.16°<length< 0.3°, 0.51° <distance< 1.1°, 
alpha< 15°). This image selection procedure ("su- 
pcrcuts") was successfully applied to repeated obser- 
vations of the Crab Nebula and has proved to be the 
most sensitive and reliable technique so far developed 
to discriminate gamma-ray events from the dominant 
hadronic background (Fegan et al. 1994). 

Motivated by the detection of Markarian 421 by 
the EGRET instrument at GeV energies (Lin et al. 
1992), observations with the 10-meter reflector were 
made initially between March 24 and June 2, 1992 for 
a total of 445 minutes on-source and an equal amount 
of time off-source during 8 clear nights. With the 
instrument operating at an effective energy threshold 
of 500 GeV, we obtained strong evidence for gamma- 
ray emission from the direction of the source (Punch 
et al. 1992). 

In the fall of 1992, several instrumental improve- 
ments were made including the recoating of most of 
the mirror facets of the 10-meter reflector. Prepara- 
tions for operation in coincidence with the nearby new 
11-meter telescope (Schubnell et al. 1992) resulted in 
changes in the trigger conditions as well. These modi- 
fications had the effect of changing the optimum selec- 
tion criteria used for selecting gamma-ray images and 
made it necessary to introduce an image size selection 
in software to adjust the energy thresholds similar to 
that used previously. This selection, which excludes 
images with a total number of less than 350 photo 
electrons, was chosen to reproduce the expected Crab 
gamma-ray rate at 500 GeV and not to optimize the 
signal. During the following winter 1992/1993, the 
AGN observation program was continued and from 
December 23, 1992 to February 26, 1993, Markarian 
421 was again observed on 10 clear nights for a total 
of 750 minutes. 

For the most recent observations between Decem- 
ber 16, 1993 and June 12, 1994, the mirror re-coating 
was completed and the trigger condition reverted to 
that of previous years. The instrument was operating 
at an energy threshold of 350 GeV, somewhat lower 
than in 1988/1989 due to the higher mirror reflec- 
tivity. Prior to the observations in May 1994, we in- 
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stalled light focussing cones in front of the PMT plane 
(Kerrick et al. 1995). This increased the number of 
Cherenkov photons collected by the photomultipliers 
and decreased the trigger energy threshold. Obser- 
vations with the light cones in place produced sta- 
ble trigger rates and first measurements on the Crab 
Nebula showed a significant increase in the observed 
gamma-ray rate. Encouraged by these results, we 
have continued to use the light cones throughout the 
entire observing season as well as in 1994/1995 with 
reproducible results on the Crab Nebula, consistent 
with a 250 GeV energy threshold. 

For the purpose of calculating the photon flux val- 
ues and to investigate possible time variations, we ref- 
erence all fluxes to the observed intensity of the Crab 
Nebula. This procedure corrects for possible system- 
atic variations which are due to the different trigger 
conditions over the entire observing period, assum- 
ing the uniformity of the Crab emission. An analysis 
of the 4-year database of Crab Nebula observations 
shows that the flux is constant to better than 10% 
(Weekes et al. 1993). This limit of variability is in- 
cluded as a systematic error of 10% to the calculated 
Markarian 421 fluxes. 

3. Results 

Because of the changes in the instrument as dis- 
cussed in the previous section we will treat each of 
the observing seasons separately. The sensitivity of 
the Whipple telescope is such that a significant flux 
of 0.3 of the Crab flux can be observed from a single 
night's observations. Therefore the dataset has been 
divided into daily intervals, with each interval span- 
ning a minimum of 1/2 hour observation. A summary 
of the observations for the three seasons reported here 
can be found in table 2. 

3.1. The 1992 Season 

The flrst observation of very high energy gamma- 
ray emission from Markarian 421 was based on data 
taken during the 1992 March- June period (Punch et 
al. 1992). We re-analyzed this dataset using the 
standard image analysis procedure but slightly modi- 
fied the algorithm to calculate the relative photomul- 
tiplier gains and excluded phototubes for which the 
pedestal distributions were outside the statistical ex- 
pectation. With the definition of the gamma-ray sig- 
nal as {Non — Noff), where Non represents the num- 
ber of selected images according to the discrimination 



procedure discussed in the previous section in the on- 
source scan and Ngff represents the number of se- 
lected images in the off-source observation, we obtain 
an average gamma-ray rate of 0.34±0.05 min~^ over 
the total 445 minutes of observation. With the statis- 
tical significance a defined as -^/A^J^T-K^Vo//, we find 
that the observed excess has a significance of 6.7a. 

Figure 1 a) shows how the difference Non — ^of / is 
distributed as a function of the orientation parameter 
alpha after applying all other selection criteria. The 
excess at small alpha values is expected for gamma- 
ray events from the source direction, while hadron- 
initiated Cherenkov events are random in orientation 
and contribute an isotropic background. Assuming 
an effective area of 3.5 x 10® cm^, the measured rate 
is equivalent to an integral gamma-ray flux of 1.59 ± 
0.20 X 10"" cm-2 sec-i above 500 GeV (see table 
3). 

Two stars (SAO 62387 and SAO 62393) with a 

combined brightness equivalent to magnitude 5 are 
located near the Markarian 421 position. This led us 
to investigate the possibility of a false signal due to 
the brightness variation in "on" and "off" fields. The 
angular separation from Mrk 421 is on the order of 
the telescope's optical resolution and thus both star 
images are always contained in the center tube. The 
most obvious test is to exclude the photomultiplier 
tube containing these stars (the center tube) from the 
analysis. No significant difference in the result was 
found. When we excluded the center tube from the 
trigger and the analysis during observations made in 
the following years, we also could find no evidence for 
a systematic effect. Additional control observations 
on two stars of similar magnitude show null results 
when subjected to the supercuts analysis. 

Figure 2 a) shows the emission on a nightly basis 
for the 1992 observing period. The errors shown are 
statistical only and do not incorporate small system- 
atic effects due to the variation in zenith angle of the 
observations. The maximum deviation (2.9 a) from 
a steady flux is found on MJD 48764 (1992 May 22 
UT). A test of the 1992 data gives a 2% probabil- 
ity that the observed flux is consistent with constant 
emission (see table 4). 

If we exclude the observation from MJD 48764 
from the sample and repeat the test on the Mar- 
karian 421 data, the probability for constant observed 
emission then increases from 2% to 40%. We empha- 
size that the observed increase in photon emission on 
MJD 48764 has a marginal significance but, in con- 
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text with the later observed day-scale flux variations, 
has to be considered as a likely flare. The apparent 
photon flux for this particular 38 minute observation 
is 5.0± 1.2 X 10-11 cm-2 sec"! above 500 GeV. 

To exclude possible systematic effects in the in- 
strument as the source for the variability in the emis- 
sion, we have analyzed contemporaneous data from 
the Crab Nebula (Figure 3a)). From this we extract 
a gamma-ray rate of 0.80 ± 0.08 min"! correspond- 
ing to an integral photon flux of 3.81± 0.03 x IQ-^^ 
cm"^ sec~i above 500 GeV. A test for constant 
flux from the Crab Nebula gives a 48% probability. 

3.2. The 1992/1993 Season 

For the 1992/1993 season, changes in the trigger 
conditions (Table 1) resulted in a loss in effective col- 
lection area in spite of the fact that a large fraction 
of the mirrors had been recoated. In order to have 
a consistent database for Markarian 421 for 1993, we 
have used only that data for which contemporaneous 
Crab observations under identical instrumental con- 
ditions are available as reference. Two instrumental 
effects had to be considered for the data analysis: 1) a 
new PMT camera with slightly different tube spacing 
was installed and 2) a persistent problem with one 
tube led to a software filtering procedure to remove 
triggers associated with that tube. 

Following the above selection, a total of 750 min- 
utes of on-source observations of Markarian 421, dis- 
tributed over 10 nights, was obtained. Applying the 
standard analysis, we obtain an average gamma-ray 
rate of 0.29 ± 0.08 min"^ with a significance level of 
3.7 cr above 400 GeV. 

From the observations of the Crab Nebula, we de- 
termined that a large portion of the trigger- induced 
background can be suppressed if we restrict the anal- 
ysis to images with a total of at least 350 photo- 
electrons. This selection increases the effective en- 
ergy threshold to 500 GeV. The corresponding on— off 
alpha distribution (with the 350 pe cut applied) is 
shown in figure lb). 

Although we have almost doubled the observing 
time compared to the 1992 season, the significance in 
the extracted signal is considerably lower due to the 
increased contamination of the selected gamma-ray 
sample with hadronic showers (Signal/Background = 
0.29 compared to Signal/Background = 0.85 for the 
1992 observations). The observed average gamma-ray 
rate is equivalent to an integral flux of 1.67 ± 0.4 x 



10 11 cm ^ sec ^ above 500 GeV, in good agreement 
with the result from the previous year if we scale the 
400 GeV 1992 result to 500 GeV, assuming an integral 
spectrum like E~^. (The assumption of an integral 
spectrum for Markarian 421 in this context is 
justified by the extrapolation of the (integral) 1992 
flux to the lower energy EGRET data points (Lin et 
al. 1994).) 

The nightly flux values were calculated to test for 
variability of the emission (Fig 2b). The maximum 
deviation {1.3a) from a steady flux is found on MJD 
49009 (1993 Jan 22 UT). A test gives a 84% prob- 
ability that the flux is constant (table 4). 

3.3. The 1994 Season 

In an effort to monitor more closely the photon 
emission from Markarian 421, extensive observations 
were taken during the period December 15, 1993 - 
June 12, 1994 (MJD 49336 - MJD 49515). The in- 
strument operated with two slightly different energy 
thresholds, depending on whether or not light fo- 
cussing cones were in place. All the scans prior to 
May 10th were taken without the light cones with 
the exception of observations on January 17 (MJD 
49369), March 5 (MJD 49415) and March 6 (MJD 
49416). For the runs without the light cones, the en- 
ergy threshold is already significantly lower {Eth=350 
GeV) than for the previous year, an improvement due 
to the higher reflectivity of the newly coated mirrors. 
The energy threshold for observations with the light 
cones is £t/i=250 GeV. 

During the night of May 15 (MJD 49487) an in- 
crease by a factor of nine in the Markarian 421 TeV 
emission was observed (Kerrick ct al. 1995a). Be- 
cause of the occurrence of this outburst we separate 
the observations into pre-burst, burst, and post-burst 
periods. 

a) pre-burst 

For observations in the standard on/off mode be- 
tween December 1993 and April 1994 we measured an 
average gamma-ray rate of 0.31 ± 0.06 min~i during 
1145 minutes, on 16 nights with a significance of 5.5 a 
(Table 2). The energy tliresliokl during this period is 
calculated from contemporaneous Crab observations 
{Rcrab = 1-89 min-i) to be 350 GeV. In addition, 
we obtained observations without the light cones in 
tracking mode on 9 nights totaling 939 minutes result- 
ing in an average gamma-ray rate of 0.17 ± 0.06. Dur- 
ing 4 nights on which light cones were in place, we col- 
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Icctcd an almost equal amount of tracking mode and 
on/off mode data (Table 2). The average observed 
gamma-ray rate is 1.59 ± 0.18. From table 2 and fig- 
ure 2c it can be seen that this high average rate is 
due to an increase in the observed emission on March 
5 and March 6. For the night of March 6, we have 
observations of the Crab Nebula which we analyzed 
to exclude the possibility of systematic effects which 
could be attributed to either the newly mounted light 
cones or the non-standard observing technique (track- 
ing) for those two nights. On March 6, 1994, we ob- 
served the Crab Nebula for 29 minutes under identi- 
cal instrumental conditions as for the Markarian 421 
observations. The detected gamma-ray rate is 2.9 ± 
1.1 min~^, in agreement with the average of 3.35 ± 
0.39 min"-'^ for all the Crab Nebula scans at the 250 
GeV energy threshold. In the absence of any other 
systematic effect (for example unstable weather con- 
ditions) the observed excess is statistically significant. 
A test applied to all pre-burst data gives a proba- 
bility of 1.3 X 10~^. The largest contribution comes 
from March 6 (see figure 2c) . The average gamma-ray 
rate for the pre-burst period (excluding March 6) for 
the combined data between December 15, 1993 (MJD 
49336) and May 9, 1994 (MJD 49481) is 0.23 ±0.04 
min~^ above 350 GeV. If we scale from the average 
of the 1992 and 1993 observing season (0.32 min~^), 
we expect a rate of 0.50 min"'^ above 350 GeV. Even 
if we allow an additional systematic error of 10%, the 
constant emission from Markarian 421 between De- 
cember 1993 and May 1994 is, by a factor of 2.2±0.5, 
significantly below that of the two previous observing 
seasons. 

b) burst 

While on May 9th (MJD 49419), Markarian 421 
was still observed at the low flux level of the previ- 
ous months, a significant increase of the photon flux 
was observed on May 15th (MJD 49487, see fig 2c), 
reaching a max;imum of 2.1x10"^° cm~^ s~^ during 
a 28 minute scan (Kerrick et al. 1995a). During this 
observation we measured a gamma-ray event rate at 
the level of 4.5±0.8 min~^, a rate comparable to the 
Crab at this energy level. The appearance of the 
Moon during the following nights made observations 
of Markarian 421 impossible again until May 30 (MJD 
49502). The TeV outburst of Markarian 421 occurred 
during a multiwavelength campaign and simultane- 
ous or contemporaneous observations are available. 
The strongest observed TeV flux precedes the obser- 
vation of an increase in the hard X-ray emission by the 



ASCA satellite by 24 hours (Takahashi et al. 1994). 
Simultaneous observations by the EGRET instrument 
did not indicate an increase in the GeV emission. A 
detailed analysis of the multiwavelength data is dis- 
cussed in a separate paper (Macomb et al. 1995). 

c) post-burst 

Following the outburst in May, we monitored Mar- 
karian 421 during the month of June when permitted 
by weather and moon. All observations were made in 
the standard on/off mode with light cones in place. 
The average observed gamma-ray rate is 0.43 ± 0.13 
min""'^ over the course of 10 nights (see figure 2c) with 
a probability of 11% for constant emission (see table 
4). 

The average integral photon fluxes measured at 250 
GeV and at 500 GeV are consistent with the extrap- 
olation of the spectrum measured between 100 MeV 
and several GeV (sec figure 4). As of the time this 
paper is written, no detection of gamma rays from 
Markarian 421 has been claimed by experiments sen- 
sitive to gamma-ray energies above 10 TeV. Reported 
upper limits on the photon emission from those ex- 
periments (Alexandreas et al. 1993, Amenomori et 
al. 1994, Catanese et al. 1994, Karle 1994, Kiihn 
1994) have been incorporated in figure 4. 

4. Discussion 

The average gamma-ray flux observed during 1992 
and 1993 corresponds to a photon luminosity of 
~ 10^^ erg s~^ above 500 GeV assuming isotropic 
emission at a distance of 124 Mpc. However, the 
fact that all identified extragalactic sources detected 
by EGRET can be associated with the blazar class 
of AGN suggests that the GeV emission is highly 
beamed, originating in the relativistic plasma out- 
flow which forms the radio jets. Beamed emis- 
sion of the gamma rays furthermore overcomes the 
strong photon-photon absorption implied by the large 
gamma-ray luminosities for the case of isotropic emis- 
sion. If we adopt the picture of beamed emission for 
the TeV component, the true VHE gamma-ray lu- 
minosity decreases by a factor of ~ 10~^ for beam 
opening angles of the order of a few degrees. 

Much attention has been devoted lately to the dis- 
cussion of possible absorption of the TeV photons in 
the intergalactic radiation field (Stecker, De Jager, & 
Salamon, 1992, Salamon, Stecker, & De Jager, 1994) 
i.e. gamma rays with energies of 0.5 TeV - 5 TeV 
interacting with soft 0.05 - 0.5 eV photons to pro- 
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duc;c' e+c^ pairs. For a relatively nearby source like 
Markarian 421 (z=0.031) this is a rather small effect. 
(The possible presence of an absorption effect in our 
original 1992 data (Mohanty et al. 1993) has been 
used by De Jager et al. (1994) to determine a value 
for the infrared energy density. A more conservative 
approach has been adopted by Biller et al. (1995) to 
establish an upper limit to the energy density.) 

The question arises whether or not this process 
could be occurring in the vicinity of the AGN, that 
is to say, can gamma rays with TeV energies escape 
the photon field close to the central engine of the 
AGN? For some models this requirement is a severe 
obstacle as explained for instance by Coppi, Kartje, 
& Konigl (1993). While their theory nicely models 
the non-thermal emission component by upscattering 
radiation (ambient or synchrotron-self- Compton pro- 
duced) by a relativistic electron/proton beam, they 
need to place the TeV emission regime further from 
the core than the GcV emission in order to avoid ab- 
sorption in the region of high energy density at small 
distances. In simultaneous observations during the 
May 1994 TeV outburst (Macomb et al. 1995), the 
EGRET instrument did not observe an increase in 
the photon emission from Markarian 421 above the 
previously detected level (Lin et al. 1994). 

Many of the various models that have been pro- 
posed to explain the origin of the gamma-ray emission 
favor production of the gamma-rays by the inverse- 
Compton process. Beamed relativistic electrons in 
the jet upscatter low energy photons from synchotron 
radiation (Ghisellini & Maraschi, 1989; Marscher 
& Bloom, 1992) or other ambient sources (Bland- 
ford 1993; Dermer, Schlickeiser, Mastichiadis, 1992; 
Sikora, Begelman, Rces, 1994) to gamma ray energies. 
These models imply a strong correlation between the 
production of radiation across the non-thermal wave- 
band. Any change in the electron flow or electron dis- 
tribution translates into a variation in the emission 
at other wavelengths. Differences in the relation of 
variability at given frequencies reflect distinct produc- 
tion regions. The inhomogcncous synchrotron-self- 
Compton model by Maraschi. Ghisellini, and Celotti 
(1992) for instance places the production of gamma- 
rays at the inner part of the jet while hard X-rays (as 
well as radio and IR emission) are produced in the 
outer jet part and thus fluctuations in the gamma- 
ray band are predicted to lead X-ray variability by 
typically a day. 

In a similar model by Sikora, Begelman, and Rees 



(1994), the entire gamma-ray spectrum is produced 
within the same jet region through Comptonization 
of ambient radiation (from dust near the jet) by 
relativistic electrons. The authors point out that 
the problem of TeV opaqueness can be overcome if 
the seed photons come from thermal IR emission by 
dust rather than UV photons, supported by the fact 
that while thermal UV excess is common for OVV 
quasars, it is only observed in very few BL Lac ob- 
jects. This hints that the non-observation of EGRET 
detected sources other than Markarian 421 (Kerrick 
et al. 1995b) above a few hundred GeV may be a com- 
bination of an intrinsic source feature (BL Lac's vs. 
OVV's for example) and the increasing IR absorption 
towards higher redshifts. Clearly, this needs more ob- 
servational evidence, some of which can be provided 
by future measurements of AGN emission spectra in 
the (10-100) GeV regime where both processes will 
potentially cut off the power-law spectrum observed 
at high energies. 

In contrast to the leptonic models is the approach 
by Mannheim (Mannheim 1993; Mannheim & Bier- 
mann 1992) in which a shock accelerated ultra-rela- 
tivistic proton population in the jet photoproduces 
tt'^'s and the observed gamma rays are then produced 
by synchrotron cascade reprocessing. In this model, 
the gamma-ray emission observed at GeV energies 
extends naturally into the TeV regime. No predic- 
tions about the correlation of variability in the TeV 
gamma-ray emission to the GeV and X-ray emission 
are made but it was pointed out by von Montigny et 
al. (1995) that a disturbance in the proton population 
can translate to large gamma-ray flares downstream 
from the cascade origin. 

Although the Whipple collaboration has recently 
reported a somewhat weaker TeV AGN source, Mar- 
karian 501 (Quinn et al. 1995), Markarian 421, is, so 
far, the best probe to test models at the highest acces- 
sible energies. The detection of a TeV component in 
the Markarian 421 emission as well as the observation 
of day scale variability already constrains models of 
gamma-ray production in relativistic jets. Additional, 
extended simultaneous observations of Markarian 421 
over the optical, UV, GeV, and TeV energy regime 
are crucial to establish the correlations of non-thermal 
photon emission with the ultimate goal to understand 
the gamma-ray production and emission process in 
blazars. 
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TABLE 1 

Instrument Parameter 



Season 


Mode 




Collection Are; 


I Trigger 






[GcV] 






1991-1992 


on/ofT 


500 


3.5x10'' 


ni=2,n=50 


1992-1993 


on /off 


400 


2.2x10* 


m=3,n=36 


1993-1994 


on/off, tracking 


350 


3.5x10* 


m=2,n=30-40 


1993-1994 (LC) 


on/off, tracking 


250 


3.5x10* 


m=2,n=30-40 



TABLE 2 



Observation Summary 



Observation Period 




Mode 


Duration 


Nights 


Ganinia-Ray Rate 




[GeV] 




[min] 




[min-i] 


Mar 24, 1992 - Jun 2, 1992 


500 


on/off 


445 


8 


0.33±0.05 


Dec 23, 1992 - Feb 26,1993 


400 


on/off 


750 


10 


0.29±0.08 


Dec 15, 1993 - May 9, 1994 


350 


on/off 


1145 


16 


0.31±0.06 


Jan 10, 1994 - May 3, 1994 


350 


tracking 


939 


9 


0.17±0.06 


Jan 17, 1994 - June 12, 1994 


250 


on/off 


665 


15 


0.88±0.14 


Mar 5, 1994 - May 30, 1994 


250 


tracking 


516 


5 


1.47±0.20 



TABLE 3 

Average Season Fluxes 



Observing Season 




Flux [cm~'^s~'] 


1991/1992 


500 


1.59±0.20 xlO" 


-ii 


1992/1993 


500 


1.47±0.30 xlO" 


-11 


1993/1994 pre 


350 


1.18±0.20 xlO^ 


-11 


1993/1994 burst 


250 


21.90±3.80 xlO 


-11 


1993/1994 post 


250 


2.48±0.60 xlO" 


-11 



TABLE 4 

Probability for Constant Emission 



Observing Season 


Markarian 421 




Crab Nebula 




Prob. 




ndf 


Prob. 


1991/1992 


1.9 x 10"^ 


16.8 


7 


0.48 


1992/1993 


0.84 


5.0 


9 


0.66 


1993/1994 pre-burst 


1.3 X 10"® 


97.4 


28 


0.78 


1993/1994 post-burst 


0.11 


15.6 


10 


no data 
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Figure 1. On-source minus ofF-source orientation angle (alpha) distributions for the three viewing periods 1992 
(a), 1992/1993 (b), and 1993/1994 (c). The net excess in gamma-ray like events accepted by the selection procedure 
is indicated by the shaded area. The dashed line shows the average net gamma-ray rate for images with large alpha 
angles. This rate is expected to be zero. A one sigma statistical error bar is always shown for the bin 85—90°. 
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Figure 2. The observed fliixes from Markarian 421 for the three viewing seasons, (a) 1992, (b) 1993, and (c) 
1994. Each data point corresponds to one night of observation. Different symbols indicate different instrument 
parameter and operation mode: solid circle - on/off mode, no liglit cones; open circle - tracking mode, no light 
cones; solid box - on/off mode, light cones; open box - tracking mode, light cones. 



12 



> 
o 



15 



10 



: 1991/1992 


1 1 1 1 1 








+ 


+ 






1,1,1 







8580 8600 8620 8640 8660 8680 8700 8720 8740 

Modified Julian Day - 40000 



20 - 



10 - 




9320 9340 9360 9380 9400 9420 

Modified JuUan Day - 40000 



Figure 3. The observed fliixes from the Crab Nebula for the viewing seasons 1992 (a), and 1994 (b). Symbols 
as in Figure 2. 
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Figure 4. Integral photon spectrum for Markarian 421. The error box for EGRET was converted from the 
differential fit to the observed data (Lin et al. 1993). Upper limits are from air shower experiments: HEGRA (Karle 
1994; Kiihn 1994); Cygnus (Alexandreas et al. 1993); Tibet (Amenomori et al. 1994); CASA-MIA (Catanese et al. 
1994). The Whipple data points arc the average integrated flux for 1992 and 1992/1993 (open circle), the pre-burst 
average flux in 1994 (open triangle), and the flux observed during the burst in May 1994 (solid circle). 
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